D° - D° mixing in Y(1S) -> D°D° decay at Super-B 
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D° — D° mixing and significant CP violation in the charm system may indicate the signature of new 
physics. In this study, we suggest that the coherent D°D° events from the decay of T(1S) — > D°D° 
can be used to measure both mixing parameters and CP violation in charm decays. The neutral 
D mesons from T(1S) decay are strongly boosted, so that it will offer the possibility to measure 
the proper-time interval, At, between the fully-reconstructed D° and D°. Both coherent and time- 
dependent information can be used to extract D° — D° mixing parameters. The sensitivity of the 
measurement should be improved at B factories or a super-B. 
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Due to the smallness of the Standard Model (SM) 
AC = amplitude, D° — D° mixing offers a unique 
opportunity to probe flavor-changing interactions which 
may be generated by new physics. The most promis- 
ing place to produce D°D° pairs with low backgrounds 
is the ^(3770) resonance just above the D°D° thresh- 
old. The current experiments, such as CLEO-c and BE- 
SIII [lj], are all symmetric D meson factories, on which 
the time information cannot be used. It is very hard to 
build an asymmetric r-Charm factory in order to sepa- 
rate the two D° decay vertices since we need a strong 
boost of the D meson to measure the mixing parame- 
ters. Although the time-dependent analyses have been 
done at B factories, the D mesons produced there are 
incoherent. In T(1S) — ► D°D° decay, both D mesons are 
strongly boosted in the rest frame of the T(1S) with the 
Lorentz boost factor of (((3j)d = 2.33), precise determi- 
nation of the proper-time interval (At) between the two 
D meson decays is available. Both coherence and time 
information are essential to measure the D° — D a mixing 
and CP violation. 

In this paper, we consider the possible observations of 
D° - D° mixing and CP violation in the T(15) -► D°D° 
decay, in which the coherent D°D° events are generated 
with strong boost. Here we assume that possible strong 
multiquark effects that involve seaquarks play no role 
in Y(IS) -> D°D° decays Q. The T(1S) decays will 
provide another opportunity to search for D Q — D° mix- 
ing and understand the source of CP violation in charm 
system. The amplitude for T(1S) decaying to D°D° is 
(D°D°\H\T(1S)), and the D°D° pair system is in a state 
with charge parity C = — 1, which can be defined as Q 



1 

3 



\D U )\D U 



(1) 



Although there is a weak current contribution in 
Y(IS) — > D°D° decay, which may not conserve charge 
parity, the D°D° pair can not be in a state with C = +1. 
The reason is that the relative orbital angular momen- 
tum of D°D° pair must be I = 1 because of angu- 



lar momentum conservation. A boson-pair with / = 1 
must be in an anti-symmetric state, the anti-symmetric 
state of particle-anti-particle pair must be in a state with 
C = -l. 

We shall analyze the time-evolution of D°D° system 
produced in T(1S) decay. 

In the assumption of CPT invariance, the weak eigen- 
states of D° - D° system are \D L ) = p\D°) + q\D°) and 
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p\D° 



q\D°) with eigenvalues /il = tul — ^ L 



and fin — fnH — 2^^' res P ec ti v ely, where the mr, and 
Tl (ttih and Th) are the mass and width of "light (L)" 
L>°J "heavy (H)" D°) meson. Following the T(1S) -> 
D°D° decay, the D° and D° will go separately and the 
proper-time evolution of the particle states l-^p^ygC*)) 

and |-Dp hys (£)) are given by 
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with definitions 
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Am = mn — mr,, 



, Ar = T H - T L 
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Note that here Am is positive by definition, while the 
sign of Ar is to be determined by experiments. 

In practice, one define the following mixing parameters 
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Then we consider a D°D° pair in Y(1S) decay with defi- 
nite charge-conjugation eigenvalue. The time-dependent 
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wave function of D°D° system with C = — 1 can be writ- 
ten as 



\D 



phy S (ki,*i))l^ hys (k 2) t 2 ))], (6) 



where ki and k2 are the three-momentum vector of the 
two D mesons. We now consider decays of these corre- 
lated system into various final states. An early study of 
correlated D°D° decays in to specific flavor final states, at 
a r-Charm factory, was carried out by Bigi and Sanda 
Xing [5j had considered time-dependent decays into cor- 
related pairs of states at -0(3770) and i/'(4140) peaks. The 
amplitude of such joint decays, one D decaying to a fi- 
nal state /i at proper time ii, and the other D to f 2 at 
proper time t 2 , is given by [5| 



A{T(1S) 



^phys^phys ~> fcfo) = 
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{a-[g-(t 1 )g + (t 2 )-g + (h)g-(t 2 )} + 
o+[ff_(ti)ff_(t 2 )-fl + (ti)5 + (*2)]} (7) 



interference between these two cases. In the following 
discussion, we define 



R(h,h;t) 



dT(T(lS) -> ^ phys ^ phys - hh) 



dt 



(12) 



The time-dependent rate expression simplifies if one of 
the states (say, f 2 ) is a CP eigenstate S v with eigenvalue 
j] = ±: 

R(h,S v ;t) = Af\A s f\A fl \ 2 e- r ^ x 
[2|A A +7 ? | 2 cosh(yrt)- 
27 7 (|A /l +7 7 | 2 )sinh( 2/ R)], (13) 

where A Sv = (S v \H\D a ), and we have used CP\D°) = 
— \D°) and As, = — 77 = T by neglecting CP violation in 
decay, D° - D° mixing and the interference of the decay 
with and without mixing. 

Now we consider the following cases for the D me- 
son decays to various final states, such as semileptonic, 
hadronic, and CP eigenstates. 

(1) (l-X+,K+7T-;t): 



where 



a+ = A h A h - A h A h = A h A h -{X h - A/J, 
a_ = V -A h A h - ^A h A h = A h A h ^{\ - X h X h ), 



with A h = (fi\H\D ), A fi = {fi\H\D°), and define 
_ q(M\D°) _ qA fi 



A/, 
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q(fi\H\D°) qA 



(10) 
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In the process e + e~ — > T(1S) — ► D°D° the central-of- 
mass energy is far above the threshold of D°D° pairs, so 
that the decay-time difference (t=At_ = (t 2 — ti)) be- 
tween -Dpjjyg - > fi an d -^phys ~~ * $ 2 can ^ e measure d eas- 
ily. From Eq. ((7|), one can derive the general expression 
for the time-dependent decay rate, in agreement with [f|: 



dT(Y(lS) -> D°, D°, hh) 

v v ; Phys phys JUZ > = u- e -T\t\ x 

dt 

[(\a+\ 2 + |a_| 2 )cosh(ylY) + (|a + | 2 - |a_ | 2 )cos(irrt) 
-2TZe(a* + a-)smh(yTt) + 2Zm(a* + a-)sm(xTt)], (11) 

where Af is a common normalization factor, in Eq. (jlip . 
terms proportional to |a+| 2 are associated with decays 
that occur without any net oscillation, while terms pro- 
portional to I a, 1 2 are associated with decays following a 

net oscillation. The other terms are associated with the 



R(rX+,K+w-;t) = Af\Ai\ 2 \A K+ , 



e- r l*l x 



((l+^A K+ ^| 2 )cosh( y rt) 
(1 - \X K +^-\ 2 )cos{xTt) 
2TZe(Xx+Tr- )smh(yTt) 
2Tm(X K + n - )sin(xTt)) . 



(14) 



(2) (l+X- ,K-TT+:t): 



R{l+X-,K-ir+;t) = U\A t \ 2 \A K -^ 



e" r l*l x 



((l + |A K - ff+ | 2 )cosh( y n) 
- (l-\\K-«+\ 2 )<xx(x£t) 
+ 2TZe(X K - 7T +)smh(yTt) 
+ 2Xm(X K - n +)sm{xTt)). 



(15) 



(3) (l+X-,K+ir-;t): 
R(l + X-,K + n-;t) 



Af\A t \ 2 \A K+7 



| 2 e- r l*l x 



((l + \X K +n-\ 2 )cosh(yTt) 
+ (1 - \X K+7T -\ 2 )cos{xrt) 
+ 2Ke(X K+w -)smh(yTt) 

- 21m(X K+7T -)sin(xTt)). (16) 

(4) [l-X+,K--K + \t): 

R{rX-,K-v+;t) = AA|A ; | 2 |A A - T+ | 2 e- r|t| x 
((l + |A K - T+ | 2 )cosh(yR) 
+ (1 - |Ax-^+| 2 )cos(afrt) 
+ 27ee(A K - 7r +)sinh(yn) 

- 21m{X K - 7I +)sm(xTt)). (17) 
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(5) (ifx^j^x^ty. 

(cosh(yR) + cos(aTt)), (18) 

where h and I2 could be electron or muon. 

(6) (l±XT, 

iZ(i = R(l + X 7 S n ;t)\qo p 

= M\Ai\ 2 \As„\ 2 e- r ^ x 



[2cosh(yTt) - 2rfRe ( - ) sinh(yn) 



(19) 



-2rfEm ( - ) sin^rt)], 



where q p indicates the exchange of q and p, and 
\q/p\ = 1 is taken. 

(7) (K-n+,S n ;t): 

For this case, it is the same as the result in Eq. (p~3|) for 
fx = K~tt + when CP violation is neglected. 

(8) {K-TT + ,K+n~;t): 

For the given final states /1/2 = (K~ir + )(K + ir~), the 
situation becomes more complicated, one can obtain the 
following expression after a lengthy calculation: 

R(K-TT + ,K + 7r-;t) = Af\A K - v +A K+7r -\ 2 e~ m x 
- A K - 7r + A^+tt- I 2 + |Ak-tt+ - X K+7r - | 2 )cosh(yrt) 
- ^K-7i+ A K+7r - 1 2 - |A A -- 7r + - A^c+tt- | 2 )cos(xrt) 
+2(TZe(X K -^+X K +^- - l)TZe(X K - 7I + - X K+7T -) + 
Im(X K - 7T +\ K+7T -)Im(X K - 7r + - A i f+ 7r -))sinh(yrt) 
-2(TZe(X K -^+X K+1T - - l)lm{X K - 7! + - X K+7r -) - 
lm(X K -„+X K +„-)ne(X K - 7r + - X K+ir -))sin(xTt)]. 

(20) 

(9) (K-TT + ,K-TT+;t): 

R{K-Tr+,K-ir+;t) = ^\A K - w+ \ i \^-\ 2 e- r ^ 

x \\ 2 K -„+ - l| 2 [cosh(yrt) - cos(xR)]. (21) 

Mixing is the necessary condition for this process to oc- 
cur. 

(10) (K+TT-,K+TT-;t): 

R(K+n- , K+7T- ; t) = N\A K+n - | 4 |2| 2 e - r l*l 

P 

x \X 2 K+7r - - l| 2 [cosh(yrt) - cos(irrt)] (22) 



(11) (l+X-,l+X-,t): 



R{l+X-,l+X-,t) = Ne- m 
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irx- 



>\A 



nx- 



x (cosh(yrt) — cos(a;rt)) . 



(23) 



(12) (l^X+,l^X+,t): 
R(l^X + ,l^X + ,t) = Afe- m 



\Arx+\ 2 \A-x+\ 2 



x (cosh(yrt) — cos(xTt)) . 



(24) 



In deriving the above formulas from Eqs. (fl4|) to (I24|) . 
we have assumed that: (1) AQ — AC rule holds, A t - = 
(l-X+\H\D°) = A l+ = (l+X-\H\D°) = 0; (2) CPT 
invariance holds. The results in Eqs.f [T4]). p5[) ([TBI) and 
(fT7|) are in agreement with those in Ref. [5J]. 

In order to simplify the above formula, we make the 
following definitions: 



2 = (1 + A M )e- i0 , 
P 



(25) 



where (3 is the weak phase in mixing and Am is a real- 
valued parameter which indicates the magnitude of CP 
violation in the mixing, and for / = K ~7r + , we define 



.4 



K-tv+ 



A 



Arc- 



ire 



K--K+ 



A 



(26) 



where r and a (r' and a') are the ratio and relative phase 
of the doubly Cabibbo-suppressed (DCS) decay rate and 
the Cabibbo-favored (CF) decay rate. Then, X K -^+ and 
Aa+tt- can be parameterized as 



Aa-x+ = ~Vr(l + A M )t 



\A+7T~ 



-i{a+0) 



-i(a'-P) 



l + A 



M 



(27) 



(28) 



In order to demonstrate the CP violation in decay, we 

define \fr = \j Rr>- -. — and \fp = \J Rd{1 + Ad)- 

1 -I- Ad 

Thus, Eqs. ([27]) and |28|) can be expressed as 



K-TT + 



\A+7 



Rd 



-l + A 



l + A, 



M e - 4 (5+0) 



'Rd 



A + A 



l + A 



d e -i(S-4,) 
M 



(29) 



(30) 



Oi + Oc 

where S = — - — is the averaged phase difference be- 



tween DCS and CF processes, and 



a — a 



■0- 



We can characterize the CP violation in the mixing 
amplitude, the decay amplitude, and the interference 
between mixing and decay, by real-valued parameters 
Am, Ad, and 4> as in Ref [7[. In the limit of CP 
conservation, Am, Ad and <f) are all zero. Am = 
means no CP violation in mixing, namely, \q/p\ = 1; 
Ad = means no CP violation in decay, for this case, 
r = r' = R D = \A K -^+ /A K - v +\ 2 = \A K+n - /A K+V - | 2 ; 
= means no CP violation in the interference between 
decay and mixing. 
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Taking into account that Xk--k+, ^-k+w- ^ 1 an d 
y <C 1, keeping terms up to order x 2 , y 2 and i?u in 
the expressions, neglecting CP violation in mixing, de- 
cay and the interference between decay with and without 
mixing (Am = 0, Ad = 0, and = 0), expanding the 
time-dependent for xt, yt < T -1 , we can write the results 
from Eqs. (H) to (JUJ) as 

(1) (l-X + ,K + TT-;t): 



21 i i2„-r|t| x 



R(l-X+, K+TT-;t) = Af\A l \ 2 \A K+ ^- \-e 
(2R D - + R M T 2 t 2 ), 



(31) 



y 



is mixing rate, and y' = ycosd 



where Rm = - 
xsin^. 

(2) 

i?(? + X-,/i-7r + ;t) =M\A l \ 2 \A K -^\ 2 e- r ^ x 
(2i? D - 2^R D ~y'Tt + R M T 2 t 2 ), (32) 

since the mixing in neutral D is tiny, it is much more 
likely that x 2 , y 2 <C Rd, cases (1) and (2) can be used 
for measuring Rd- 
(3) {l + X-,K+it-;t): 



R{l+X-,K+w-;t) = XT\A l \ 2 \A K+7r - \ 2 e 
(2 - 2^/R^(ycos(S) + xsm(5))Tt 



-r\t\ 



y 



-T t 



(33) 



(4) (l-X+,K--K+\t): 

R(rX-,K-ir+;t) = Jsf\Ai\ 2 \A K -^ \ 2 e~ r ^ 
(2 - 2^/RD~(ycos{5) + xsin(<5))rf 

+V— —T 2 t 2 ) = R(l+X-,K+7r-;t). 



(34) 



In the limit of no CP violation, case (3) is the same as 
(4). 

(5) (lfX^,qX ± ;t): 



RilfX^J^X*; t) = M\A h \ 2 \A l2 | 2 e- r l*l 



(2 



(6) (!±,S v ;t): 



-Th 2 ) 



(35) 



R(l ± ,S r ,;t)=M\A l \ 2 \A Sr} \ 2 e- r m x 

(2 - 2r](ycosf3 T xsin(3)Tt + y 2 T 2 t 2 ), (36) 

where y may be determined because phase (3 = 
arg[(V us V*J(V cs V u * s )} ~ 0. 
(7) (K-TT+,S v ;t): 

R(K-n+ , S n ;t) = Af\A K ^ + \ 2 \A s J 2 e- r M x 

(r) - /R^cos(5) 2 (l - rjyTt + ly 2 (Tt) 2 ), (37) 





(Rd, x, y) 




a : (0.0026,0.06,-0.02) 




b : (0.0026,0.02,-0.02) 


\\ \ c 


c : (0.0026,0.06,-0.04) 


// \ \ 


d: (0.0,0.06,-0.02) 


d 
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FIG. 1: Illustrative plot of the sensitivity of the decay rate 
R(l^X^, K^^; t) to the mixing parameters of x, y and the 
ratio Rd, where 5 = 10° is taken. 



where cosS can be measured in this case by combing r\ 
— 1 and +1 final states. 
(8) (K-Tr + ,K + 7r-;t): 

R{K-ir + ,K + ir-;t) = Af\ A K - w +A K+n - \ 2 e' m x 



2 2 

(2 + V ~~ X T 2 t 2 - AR D cos(26)). 



(9) (K--K+,K--K + \t): 

R(K-tt+ ,K-7r+;t) = AflAK-^^e-^ x 
(l-2fl D cos(2£)) x2+y2 r 2 f 2 . 



(38) 



(39) 



This process is proportional to mixing rate Rm , and can 
be used to measure the mixing parameter directly. 

(10) (K+ir-,K+Tr-;t): 

The result is the same as R(K~tt + , K~it + ; t) when CP 
violation is neglected. 

(11) (l+X-,l+X-,t): 

R(l+X-,l+X-,t) = Me- r ^\A lt x-\ 2 \Alx-\ 2 x 

^±^T 2 t 2 . (40) 

One can also definitely measure the mixing rate in the 
like-sign processes as in case (9) and (10). 

(12) (l7X+,qX+,t): 

The result is the same as (l^X^, l^X~, t) when CP vi- 
olation in mixing and decay is neglected. 

Note that in all the above cases, when CP violation in 
decay is neglected, there is = lA^+^-l. 

The experimental data from CLEO-c yield that the al- 
lowed values for the mixing parameters x and y are: y = 
-0.058±0.066, x < 0.094 @]. The ratio of the DCS decay 
rate to the CF decay rate is R D ~ {V cs V cd ) 2 » 0.0026. 
The illustrative plot of the decay rate i?(^ ± X =F , if T 7r ± ; t) 
is shown in FigQ] by taking 5 = 10°. The decay rate 
R(l^X^ , K^ir ; t) is very sensitive to the mixing param- 
eters x, y and the ratio of Rd in the region Tt ~ 1 — 6. 
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The other decay rates in Eqs. (|33[) - (|38|) are not sensitive 
to the mixing parameters and Rd, because in these de- 
cay rates, the x, y and Rp are only small correction to 
the dominant contributions. 

In the Cabibbo-Kobayashi-Maskawa (CKM) frame- 
work CP violation in the neutral D system is very small 
and can be safely neglected. However extension of the 
SM could induce new physics of CP violation [4, 8]. The 
most likely sizable effect is a possible new CP violation 
phase, 4> — a.vg(qA/pA), occurring in the interference be- 
tween mixing and decay amplitudes. Thus, in the pres- 
ence of CP violation in the interference, we can construct 
the following CP observable based on the previous calcu- 
lations. We can look at the difference between the cases 
(3) and (4), and define: 

+ _ R(l+X- } K+n-) - R(1-X+,K-k+) 
CpW ~ R(l+X-,K+n-) + R(l-X+,K-n+) ' 1 ' 

As discussed in Ref [5|, the signal is due to the interplay 
of DCS decay and mixing. With help of Eqs. (|16jl and 
(jTTJ) , we obtain 



A+- P (t) 



>(ysin<5 — xcos<5)sin</> x Ft. (42) 



Here we keep both x and y terms since the current ex- 
periments indicate that they may be at the same order, 
this is different from Ref [5J. The above asymmetry de- 
pends on non-vanishing phase 0, and also the mixing 
parameters. Within the SM, it is of order O(10~ 3 ) [§], 
which makes such an asymmetry unmeasurably tiny un- 
less there is new physics • By looking at the difference 
between R(l + , S^; t) and R(l + ,S ri ;t) in case (6), we get 
the following asymmetry 



A s v ,,s _ R(l + , S v ) - R(l ,S V ) 
cp[) ~ R(l+,S v ) + R(l-,S v y 

With the help of Eq. (JTjJJ), we obtain: 
A s £ p (t) = rjxsiup x Ft, 



(43) 



(44) 



where (3 is defined in Eq. (|25[) . This CP term depends 
on the mixing parameter x and the phase in the mixing 
amplitude. 

In experiments, the KEK-B can move to T(1S) peak 
without losing luminosity, more than 330/6 -1 data per 
year could be taken The measured cross section at 
Y(IS) peak is a R = 21.5 ± 0.2 ± 0.8 nb by CLEO 0. 
One can estimate the total T(IS) events with one year 
running of KEK-B are about 7.1 x 10 9 . While, at 
super-KEKB, about 10 11 T(1S) events can be obtained 
with one year data-taking if the design luminosity is 
8 x 10 35 [13J. More recently, a super-B factory based on 
the concepts of the Linear Collider (LC) is proposed [3], 
the low energy beam and high energy beam are 4.0 and 
7.0 GeV, respectively. The machine can run at both 



T(4S) and Y(lS) peaks with luminosity about 1.0 x 10 36 . 
Then about 10 12 Y(1S) events can be collected with one 
year's run. In Ref. [15], one had estimated the ratio of 
BK{Y{1S) -> D+D-) and BK(Y(IS) -» K+K~) as 



BK{Y(1S) -* D+D- 
BK{Y{1S) -»• K+K- 



1 

02 



» 1.0 



(45) 



where the current upper limit of B1Z(Y(1S) — » K + K~) 
is 5.0 x 10~ 4 at 90% confidence level. One can expect 
that the decay rate of T(IS) — > K + K~) is at the or- 
der of 10 -6 Qjl, so that we can estimate B1Z(Y(1S) — > 
D + D~) — 10~ 4 - 10~ 5 . At the supcr-B factory, around 
10 7 — 10 s D°D° pairs can be collected in one year's data 
taking, which is comparable with that at the BESIII with 
four years integrated luminosity (l6j 

It is known that one has to fit the proper-time distri- 
bution in experiments to extract both the mixing and 
the CP parameters, we discuss the following two cases: 
(1) Case-I: at a symmetric T(15) factory, namely, the 
Y(IS) is at rest in the Central-Mass (CM) frame. Then, 
the proper-time interval between the two D mesons is 
calculated as 



At = (r D o - r no ) 



cIPI 



(46) 



where rjjo and r^ are D° and D° decay length, respec- 
tively, and P is the three-momentum vector of D°. Since 
the momentum can be calculated with Y(IS) decay in 
the CM frame, all the joint D°D° decays in this paper 
can be used to study D° — D mixing and the CP viola- 
tion in the symmetric Y(1S) factory. (2) Case-II: while, 
at an asymmetric Y(1S) factory, the Y(1S) will be pro- 
duced with a boost. In this case, the momentum of D° 
and D° will be different each other, and one has to fully 
reconstruct at least one of the two D mesons, since the 
the proper-time interval between the two D mesons is 
calculated as 



At = r D a 



m D 



cP 



DO | 



CP 



(47) 



D°l 



where |Pd°| an d |Pgo| are momentum of D° and D°, 
respectively. In this case, the joint decays to di-lepton in 
Eq. ([TH]) will not work, since both the D mesons cannot 
be fully reconstructed. One cannot obtain the proper- 
time from such kind of experiment. 

The average decay length of the D° meson in the rest 
frame of Y(1S) is ct d o x (/3j)o° ~ 290/im. At B factories, 
such as Belle detector, the impact parameter resolution 
of vertex detector, which are directly related to decay 
vertex resolution of D°, are described in Ref. [l7], from 
which we can get that the resolution for the reconstructed 
D° decay length should be less than 100/im within the 
coverage of the detector. This means Belle detector is 
good enough to separate the two D° decay vertices, so 
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that the mixing parameters can be measured by using 
time information. 

All the results in this paper can also be applied to the 
following processes, 



r(lS) -» ir°(D Q D°) c = 



(48) 



where the D°D° pair are in a C = — 1 state, for example, 
Y(IS) -> D*°D° -> ir a D a D a . 

In conclusion, we suggest that the coherent D a D a 
events from the decay of T(1S) — > D°D° can be used 
to measure both the mixing parameters and the CP vi- 
olation in charm decays. The neutral D mesons from 
Y(1S) decay are strongly boosted, so that it will offer 
the possibility to measure the proper-time interval, At, 
between the fully- reconstructed D° and D°. Both co- 
herent and time-dependent information can be used to 
extract D° - D° mixing parameters, in which the sensi- 
tivity of the measurements could be improved by com- 
paring to the future measurement at the BESIII with the 
same amount of D°D Q pairs. 
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